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Fast frictionless dynamics as a toolbox for low-dimensional Bose- 
Einstein condensates 
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Abstract. - A method is proposed to implement a fast frictionless dynamics in a low-dimensional 
Bose-Einstein condensate by engineering the time-dependence of the transverse confining potential 
in a highly anisotropic trap. The method exploits the inversion of the dynamical self-similar scaling 
law in the radial degrees of freedom. We discuss the application of the method to preserve short- 
range correlations in time of flight experiments, the implementation of nearly-sudden quenches 
of non-linear interactions, and its power to assist self-similar dynamics in quasi-one dimensional 
condensates. 
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Counterintuitive as they are, implementations of fast 
frictionless dynamics (FFD) providing a shortcut to adi- 
abaticity in quantum systems have recently been intro- 
duced theoretically [1—3] and demonstrated in the labora- 
tory [4,5]. 

In this Letter, FFD is exploited as a tool-box to ma- 
nipulate and control low dimensional quantum gases. As 
a primary goal, we propose a method to tune the non- 
linearity of the effective low dimensional (ID and 2D) dy- 
namics of an anisotropic Bosc-Einstcin condensate (BEC) 
by engineering the time-modulation of the transverse con- 
finement. The possibility of attaining this goal by means 
of a multi-scale expansion method was recognised by Stal- 
iunas et al. for slow periodic modulations [G] and led to 
the observation of Faraday patterns in cigar-shaped BECs 
[7]- 

Our method, not restricted by adiabaticity in the trans- 
verse dynamics, provides an alternative way to implement 
a variety of schemes for the generation, stabilisation and 
control of solitons, and the study of related non-linear 
matter- wave phenomena in BEC, where it is often nec- 
essary to implement a time-dependent coupling constant 
[8] . Other situation where it can be applied is in the con- 
trolled expansion of BEC clouds [9] as in the simulation 
of cosmological analogues [10]. Moreover, we shall discuss 
how FFD can be exploited to implement nearly-sudden 
quenches of the mean-field non-linear interactions, pre- 



serve quantum correlations in time-of-flight and induce 
a self-similar expansion in interacting quasi-lD atomic 
cloud. 

FFD allows to evolve from an initial state to a final one 
without exciting the system while doing it in a given time 
t much smaller than that required for an adiabatic dynam- 
ics. FFD relies on the inversion of dynamical scaling laws, 
which can often be exploited to described harmonically 
trapped ultracold gases, such as the Calogero-Suthcrland 
model [11], Tonks- Girardeau [12,13], strongly interact- 
ing mixtures [14], Lieb-Liniger gases [15], Bose-Einstein 
condensates (BEC) [2, 16, 17], including dipolar interac- 
tions [18], and more general many-body quantum systems 
[19]. In a harmonic trap with time-dependent frequency, 
the self-similar evolution of the single-particle states <f> n 
(n = 0, 1, 2, ... ) follows from the well-known scaling law, 



tj> n {x,t) = <j> n (x/b(t),0)e tmx2 ^ 2bh - lE ^/ h /, 



(1) 



where the scaling factor 6 = b(t) is the solution of the 
Ermakov differential equation 



b + Lo 2 (t)b = Lu(0) 2 /b 3 



(2) 



satisfying the boundary conditions 6(0) = 1 and b(0) = 0, 
with E„ = hu(Q)(n+l/2), and S(t) = j*dt'/b 2 (t') [20]. 
In particular, the probability density reads \<j) n (x,t)\ 2 = 
\<j> n {x/b{t), 0)\ 2 /b(t). The essence of FFD is to exploit the 
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existence of the self-similar dynamics to force a desired tra- 
jectory b(t) and invert Eq. (2) to determine the required 
modulation of the control parameter uj(t) [1,2]. 

Effective one-dimensional time-dependent 
Gross-Pitaevskii equation. — Our aim is to find a 
ID effective non-linear Schrodinger equation or Gross- 
Pitaevskii (GPE) equation for a BEC in a elongated trap 
in which the transverse confinement is modulated in time, 
V cx (x,t) = V cx 0) + f [to x (t) 2 x 2 + Lo y (t) 2 y 2 }. We will do 
so exploiting the scaling law in Eq. (1). We start with 
the 3D timc-dcpcndcnt GPE which governs the dynamics 
of the order parameter '5 = ^(x, t) 

8<S> fr 2 

ih ~dt = [ - 2^ A x + V cx (x,t) + . 93 d|*| 2 ]*, (3) 

with the normalisation condition J dx|\&(x, t)\ 2 = 1, and 
where 03D = 4nh Na , N is the number of atoms in the 

condensate of mass to, and a the s-wave scattering length 2 ^ L I ^» J In I ^ j 
which determines the healing length £ = 1/ s/%wn\a\. It is 
convenient to introduce the characteristic length in each 
direction, a,j = ^/hj '(toWj(O)) for j = x, y, z, in terms of 
which the density scales as n ~ N/ Ylj a j ■ For tight trans- 



warrants a smooth modulation of the transverse trapping 
frequency and avoids abrupt changes at t = 0. Nonethe- 
less, abrupt changes of the control parameter are often 
exploited in bang-bang control methods [ ] . 

As pointed out by Olshanii, the transverse confinement 
can severely modify the scattering properties, and in par- 
ticular, lead to a confinement-induced resonance of the 
form giu — > .9113(1 — Ca/a^)^ 1 , with C = 1.4603... 
[25,26]. Eq. (5) holds far away from confinement in- 
duced resonances, i.e. \a\ <g; {a x ,a y }. In addition, low- 
dimensional BECs generally exhibit phase fluctuations 
which can lead to deviations from the mean- field [27, 28] . 
Even at zero-temperature, quantum fluctuations suppress 
the off-diagonal long-range order, in the sense that the 
reduced one-body density matrix pi(x,x') — > x( x )*x( x ') 
with x{ x ) = y/n{x) exp(— To) for \x — x'\ ~ R being R the 
size of the system and n(x) the local density. The quan- 
tum suppresion parameter for D = 1 is given by T\ = 



where ui\ 



verse confinement {uj x 



UJ, 



3> uj z ), the kinetic energy 



in the transverse direction governs over 2-body collisions. 
Whenever the dimensionality parameter e 2 = a x a y /^ 2 ~ 
N\a\/a z <C 1, the transverse excitations are frozen and a 
dimensional reduction of the 3D GPE is possible [21-23]. 
The 3D GPE becomes a linear Schrodinger equation for 
the radial degrees of freedom. We next use the ansatz 
Vf^x, t) = 3>o(x,y,t)i[j(z,t) and consider the ground state 
condensate wave function in the transverse direction at 
t = to be that of an unperturbed 2D harmonic os- 
cillator, $0(^,2/, t = 0). It follows from Eq. (1) that 
|$o(x,y,i)| 2 = \MxMt),y/b y (t),t - 0)\ 2 /(b x (t)b y (t)). 
After dimensional reduction, the following effective ID 
GPE with a time-dependent non-linearity is derived, 

lh ^T^ = [Hz+gMm(z,t)\ 2 ]ij(z,t), (4) 

where H z = -^J^r + V ox (z) — a(t), and the term a(t) = 

can be removed with a unitary transformation. In the 
following, we focus our attention of the effective coupling 
constant defined as 



W 2 (0)+ W 2 (0) [29]. 
The applications discussed in the rest of the manuscript, 
involve the monotonic expansion of the cloud in the trans- 
verse direction (b(t > 0) > 0) which diminishes the role of 
phase fluctuations. It will suffice for us to focus on situa- 
tions where phase fluctuations are negligible in the initial 
state so as to prevent the formation of density ripples for 
t > [30], i.e. Ti <C 1. However, should one be interest in 
increasing the effective interactions to a maximum value 
<7id(t), the condition for phase fluctuations to be negligi- 
ble becomes Pi(7id(t)/<7id(0) <C 1. We note as well that 
the scaling function b(t) > 0, so that it is not possible to 
change the sign of the interactions as with Fcschbach or 
confinement-induced resonances [26]. This technique, is 
therefore restricted to tune the amplitude of the coupling 
constant. 

Let us now consider an isotropic transverse confine- 
ment such that uj±(t) = ui x {t) = io y {t) (b — b x = b y ), 
and assume we are interested in a given time dependence 
9ir>(t) = 9iD(0)/b(t) 2 . This can be engineered by a trap 
modulation of the transverse trapping frequency given by 



= w±(o) 2 



gig 
fflD(0) 



2 / \ 2 

, 1 9id 3 / .giD 



2 gin 



9m 



(7) 



9m(t)=g 3 n dxdy\$ (x,y,t)\ l 



2h 2 Na 



satisfying the conditions 
2h 2 Na 



3id(0) 



ma x a y 



, . 9 id(0)=0, £id(0) = 



(5) 



(G) 



The third equality is optional, it follows from imposing 
the first two and the Ermakov equation with the continu- 
ity restriction uj r=x y {t = + ) = u r=x y {t = 0~). Hence, it 



where g^n = gm(t)- Note that this expression is not 
positive-definite, and Wj_(t) might involve imaginary fre- 
quencies as discussed below. It is straightforward to verify 
that for ballistic expansion along the transverse degrees 
of freedom, which leads to a polynomially decaying non- 
linearity gm(t) = ,9id(0)/[1 +w]_(0)£ 2 ] [31], the required 
transverse frequency is indeed u>\(t > 0) = 0. Moreover, 
the existence of the invariant of motion [1] associated with 
the Ermakov equation underlying the derivation of Eq. (7) 
allows to drive the transverse frequency without fulfilling 
the adiabaticity condition uj±(t)/uj 2 _ (t) <C 1 as long as the 
decoupling of the transverse degrees of freedom hold. This 
allows ultimately to engineer fast modulations of the ef- 
fective low-dimensional coupling constant gm(t) even in a 
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time scale comparable to wj 1 . Note nonetheless that the 
excitation of parametric resonances can distort the dy- 
namics [6], as in the experiment [7] were the observation 
of Faraday waves was reported. 

Preserving short-range correlations in time-of- 
fiight measurements. — Time of-flight expansions con- 
stitute an essential tool to study quantum correlations in 
ultracold gases. Under the assumption of ballistic dynam- 
ics it is possible to relate the asymptotic density profile 
with the momentum distribution of the initial (trapped) 
state. The slow power-law decay of the interactions 
9iT>(t) = <7id(0)/[1 + w 2 (0)i 2 ], which follows from sud- 
denly switching off the transverse potential, becomes neg- 
ligible in the time scale t ^> wj . As a result, mean-field 
interactions are non-negligible in an early stage of the ex- 
pansion, and blur short-range correlations in a length scale 
Sx ~ c/loj_ where c is the speed of sound. To avoid this, 
one would like to suppress the role of interactions in a 
faster time-scale r. To this aim, consider the reverse engi- 
neering of the time-modulation of the transverse trapping 
frequency that lead to 

9iu(t) = 5iD(0)sech-, (8) 

T 

which satisfies the first two boundary conditions in Eq. 
(6). Using Eq. (7) we find that it can be induced by 
wi(t) = (W(0) 2 - J-^sech 2 ^ - 3 i 7 . For short decay 
times r, this trajectory leads to negative values of U)\{t) 
associated with purely imaginary frequencies. The physi- 
cal implementation of uij_(t) < requires that the trans- 
verse confining potential becomes an expulsive barrier, 
pushing the atoms away from the longitudinal axis of the 
cloud. This is a general feature of trajectories given by 
Eq. (7) in combination with the first two boundary condi- 
tions in Eq. (6), and is often required for modulations of 
the coupling constant gin and transverse density in a time 
scale smaller than cjJ . For the time dependence in Eq. 
(8), whenever r < y= 1 , ui±(t) 2 < for all * > 0. For 

V 2u± (0) 

T > ^ 1 (o) ' L0 ^) 2 > at early stages and becomes neg- 
ative only after t* = rarcsech ( . 1 I . Expulsive 

y ' VV4"-l(0) 2 t=-1 J ^ 

potentials have already been used in the laboratory, for 
instance in the study of bright solitons [32], and can be 
generated in a variety of ways as discussed in [1], More- 
over, provided that <7id(0) ^ 0, the implementation of the 
trajectory generally requires sudden jumps of w±(t) as in 
bang-bang control methods [24]. 

Inducing FFD in quasi-lD Bose-Einstein con- 
densates. — The existence of scalings laws leading to 
a self-similar dynamics has motivated FFD proposals for 
supcrfast expansions providing a shortcut to adiabatic- 
ity [1,2], as recently demonstrated in the laboratory both 
with ultracold gases [4] and BECs [5]. When such tech- 
niques are extended to quasi-lD BEC, the implementation 
of a self-similar dynamics generally requires to modulate 




t/T 

Fig. 1: Shortcut to adiabaticity for a quasi-lD interacting BEC. 
Outside the Thomas-Fermi regime, the required self-similar 
dynamics in the axial direction is assisted by a simultane- 
ous modulation of the transverse confinement. The parame- 
ters are u x {t) = 0.01w z (0) with b(r) = 10, u±(0) = 20u z (0), 
t = 4oj 3 (0) _1 . The inset shows the smooth evolution of the 
corresponding scaling factors. 

in time both the axial trapping frequency and the coupling 
constant. As a result, proposals for FFD based on scal- 
ing laws become particularly challenging. In the following, 
it is shown that a modulation of the transverse trapping 
frequency can be used to tunc the axial effective coupling 
constant without the need to use a Feschbach resonance. 

For the ID GPE with a time-dependent axial harmonic 
trap 

, « r h 2 d 2 muj z (t) 2 z 2 . . t .fl"i 

it is possible to exploit the scaling law for ^> = ty(z,t) 
in Eq. (1) for the order parameter, with the chemical 
potential /i playing the role of the eigen-energy [2]. The 
condensate wavefunction evolves self-similarly according 
to 

only if the axial scaling factor & z , obeys the Ermakov equa- 
tion b z +uj 2 z (t)b z = uj 2 z (0)/b 3 z (and b z (0) = 1, b z (Q) = 
0), and the time-dependent non-linearity of the form 
9io(t) = 9id(0) /b z (t) is implemented. The required time- 
dependence of giD(t) could be achieved cxploting a Fes- 
chbach resonance to tune g?,D- Nonetheless, for a narrow 
resonance a fast change of the coupling constant -requiring 
a high control of the external magnetic field- can be exper- 
imentally challenging [26]. Alternatively, under isotropic 
transverse confinement, an axial self-similar dynamics un- 
der u) z (t), can be assisted by keeping g3D constant and 
changing the transverse trapping frequency along the tra- 
jectory 
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which implements the required time-dependent non-linear 
coupling. This trajectory induces a frictionless dynamics 
in the transverse direction (essentially a free harmonic os- 
cilator), which modulates the three dimensional density, 
and ultimately the effective axial non-linearity in the re- 
quired way gio(t) = <?id(0)/& z (*) f° r the axial dynamics 
to be self-similar. Once the axial dynamics is self-similar, 
one can engineer b z (t) to drive an axial FFD. As a result, 
to perform an axial FFD of a quasi-lD BEC in a time 
r, one can proceed in the following way: i) choose the 
desired initial and final states (b z (0) = 1, b z (r)) and de- 
termine b(t) as in [1,2], ii) find the required axial trapping 
frequency 

(u) 

iii) use Eq. (10) to derive the required transverse trapping 
frequency. Implementing both Eqs. (10) and (11) leads 
to the self-similar dynamics in Eq. (9) along the designed 
FFD trajectory b z (t) . Figure 1 shows an instance of these 
trajectories. Note that thanks to the relation between the 
axial and transverse scaling factors b± (t) = b z (t) 2 , the 
self-similar dynamics along the z-axis can be assisted by a 
trajectory oj\(i) involving only real frequencies, without 
the need to implement a transverse expelling potential but 
at exceedingly short expansion times r or large expansion 
factors 6(r). 

Nearly sudden quenches. — Another useful tool to 
probe ultracold gases is the use of a sudden quench of the 
interactions. This is the case for the applications men- 
tioned above, as well as for the generation of shock waves 
and solitons, studies of relaxation dynamics and many 
other examples. Let us consider a finite time quench be- 
tween and initial 3id(0) = gi and final value giDi^) = gf 
of the coupling constant, which approaches the sudden 
limit as t — > 0. The situation resembles that of FFD. 
We require vanishing first and second order derivatives of 
gio(t) both at t = and t = r to avoid transverse excita- 
tions at the end of quench. It is convenient to consider a 
polynomial ansatz gm{t) = X)f=o a; ^' wnose coefficients 
{a{\ are completely determined by the boundary condi- 
tions, leading to the explicit form of the quench 

gi D (t) = Sl + [3s(2s - 5) + 10] 7 s 3 (12) 

where 7 = gf — gi and s = tjr. The computation of 
the corresponding trajectory u>±(t) is straightforward us- 
ing Eqs. (7) and (12) but rather lengthy as to be displayed 
here. Fig. 2 shows the solution ui±(t) to implement dif- 
ferent type of quenches. The required time to achieve a 
given ratio a = gf/gi < 1 under free evolution is tq = 
wj 1 (0)'\/7~ 1 — 1, while by implementing u±(t) it is pos- 
sible to sped up the decay by several orders of magnitude 
with a moderate modulation of the transverse confine- 
ment. Note that a = gf/gi = b±(r)~ 2 = uj±(0)/lu±(t), 
and that the amplitude of the required frequency scales 




t/x t/x 



Fig. 2: Engineering a quench of the interactions, a) Decay 
of the effective non-linear coupling constant gio(t) under free 
expansion in the transverse direction (co±(t > 0) = 0) and 
different modulations of ujj_(t) engineered according to Eqs. 
(7) and (12). b) Associated trajectories of uj±(t) for different 
quenches of g\u{t)- 

as 1/[u;_l(0)t] 2 , providing an effective lower bound to the 
achievable values of r. 

Quasi-2D condensates. — So far we have implic- 
itly focused on a cigar-shape condensate. One can simi- 
larly modulate the nonlinearity in a pancake-shaped con- 
densate under strong enough radial confinement such 
that the dynamics along the most tightly confined di- 
rection decouples from that in the BEC plane. Assume 
an oblate 3D harmonic trap with uj x 3> tj y = uj z = 
u! r . Letting ^(x, t) — <f>o{x, t)tp(y, z, t) and 4>o{x,t = 
0) = exp(— x 2 /(2a^)/('K 1 / 4 y/a^), using the scaling law for 
4>o(x, t), upon integration of the transverse coordinate, the 
effective coupling constant for a 2D cloud undergoing a 
modulation of the transverse confinement is 

92D(t) - 5 3D J dx\M^t)\ 4 - (13) 

where 3 2 d(0) = gso/ '(V^a z ), with # 2D (f = 0) = g- 2 u(0), 
<?2d(0) = 0, and as in Eq. (6), ^2d(0) = prevents discon- 
tinuous jumps of ui r (t) at t = from happening. A given 
time dependence of <?2D = <?2d(£) follows from a trajectory 

^) 2 =^(0) 2 (^) 4 + ^-2(^) 2 , (14) 

V52D(0)/ 2D V.92D/ 

which may imply imaginary frequencies associated with an 
expulsive potential, as it happens in Eq. (7) . In partic- 
ular, note that suddenly switching off the transverse po- 
tential leads to an essentially linear-in-time decay of the 
interactions 

J72D(t)=ff2D(0)/Vl + wl(0)* 3 . 

For a <72d(0 as in Eq. (8), decaying in a time scale 
r < a;" 1 , according to Eq. (14) it is found that the re- 
quired time dependent trajectory of the control parameter 
co 2 (t) = u^(0)sech 4 (£/T) — 1/r 2 < requires an expul- 
sive potential for all t > 0. We close noticing that in a 
pancake-shaped cloud, quantum fluctuations are negligi- 
ble whenever T2 = —7=— <C 1 [291 and that, should one 
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be interested in tuning the interactions to a larger value 
<?2d(t), the condition becomes r 2 g2D(T)/ff2D(0) -C 1. 

Discussion and conclusions. — In combination with 
a spatial dependence of the transverse confining poten- 
tial, FFD paves the way to control the effective coupling 
constant both in time and space. Moreover, further ap- 
plications can be envisaged such as the preparation of 
atomic Fock states by many-body atom culling methods 
[33], where starting with a large trapped atomic cloud a 
controlled increase of the interactions can provide the ex- 
pelling mechanism for the excess of atoms. 

In conclusion, we have presented a scheme to implement 
a fast frictionless dynamics of a low-dimensional Bose- 
Einstein condensate, in which spurious excitations are 
avoided without the need to fulfil adiabaticity constraints. 
Exploiting the self-similar dynamics in the strongly con- 
fined degrees of freedom, we have shown that this can 
be achieved by engineering the modulation of the trans- 
verse confinement of the cloud in an elongated trap. As 
a result, it is possible to tune the amplitude of non-linear 
interactions in these systems. We have further applied 
the method to preserve short-range correlations in time- 
of-flight, assist shortcuts to adiabatic expansions in quasi 
ID interacting BEC, and implement nearly sudden inter- 
action quenches. More generally, we argue that inverting 
the equations associated with self-similar scaling laws, al- 
lows to determine the trajectory of the control parameter 
for different processes, and constitute a powerful toolbox 
for the manipulation of ultracold atoms. 

* * * 
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